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ABSTRACT 

The amount of dust estimated from infrared to sub-millimetre (submm) observations strongly 
depends on assumptions of different grain sizes, compositions and optical properties. Here we 
use a simple model of thermal emission from cold silicate/carbon dust at a range of dust grain 
temperatures and fit the spectral energy distribution (SED) of the Crab Nebula as a test. This 
can lower the derived dust mass for the Crab by ~50% and 30-40% for astronomical silicates 
and amorphous carbon grains compared to recently published values (0.25Mq ^ 0.12Mq 
and 0.12Mq ^ 0.072Mq, respectively), but the implied dust mass can also increase by as 
much as almost a factor of six (0.25Mq 1.14M0 and O.12M0 O.IIMq) depending 

on assumptions regarding the sizes/temperatures of the coldest grains. The latter values are 
clearly unrealistic due to the expected metal budget, though. Furthermore, we show by a 
simple numerical experiment that if a cold-dust component does have a grain-temperature 
distribution, it is almost unavoidable that a two-temperature fit will yield an incorrect dust 
mass estimate. But we conclude that grain temperatures is not a greater uncertainty than the 
often poorly constrained emissivities (i.e., material properties) of cosmic dust, although there 
is clearly a need for improved dust emission models. The greatest complication associated 
with deriving dust masses still arises in the uncertainty in the dust composition. 

Key words: Stars: AGB and post-AGB, supernovae: general, individual: Crab Nebula; ISM: 
dust, extinction 


1 INTRODUCTION 

Observations suggest Lytt systems, quasars and gamma-ray burst 
hosts at high-z contain very large amounts of dust (see, e.g . 
Bertoldi et al.l 1200^ : iBeelen et al.ll200^ : iMichalowski et al1l200^ 
2010al f9). which forces models of stellar dust production and galac- 
ti c dust evolution to extremes in order to reproduce th ese results 
( towek et al.l 1 20071 : lOalL Andersen & Hiortij 1201 lal lbl iMattssoij 
IMIh . Assuming stars as the primary source of dust suggests 
the dust is due to supernovae (SNe) and a top-heavy initial mass 
function favouring the formation of SNe. This scenario can in 
principle explain the observations, provided the dust-destruction 
rate in the interstellar medium is much lo wer than expected 
based on the kinetics of SN ex plos ions (see lilDwek et alJl2()0^ : 
lOall. Andersen & Hiorthl[201 lal la: ?: lMattssonll201 1 ). 

Recent observations in the FIR/sub-mm of nearby core¬ 
collapse supernova remnants (SNRs), with ages > 2 5 years, sug¬ 
gest high dust-formation efficiencies in the SN ej ecta (iBarlow et al.l 
I 2 OI 0 I : iMatsuura et alJl^lll : iGomez et al.l 12013 . henceforth G12) 


* E-mail: larsmat@kth.se 


though these are limited in sample size. The dust masses esti¬ 
mated using canonical single or two-temperature component fits 
to the spectral energy distributions in the FIR/-submm imply dust 
masses which are uncomfortably close to (or indeed exceed) the 
am ount of metals predicte d to be ejected in these supernovae (see 
e.g. IMatsuura et al]l201 ll . G12). The Inferred dust mass depends 
on the assumed dust properties, which can thus change whether 
SNRs appears to contain unrealistically large amounts of dust or 
not. If the grains are silicates, the grain masses are high relative 
to the emissivity. For carbonaceous dust the mass/emissivity-ratio 
is much smaller, which is also the case for ice-coated grains and 
various other types of coagulates. 


Variation in dust temperatures, grain sizes and material prop¬ 
erties makes dust-mass estimates from SEDs uncertain, which 
is recognised in the literature. The fact that SEDs arising from 
dust emission often suggest there is a range of grain temper¬ 
atures in a dust population, has inspired several publications 
on how to solve the inverse problem of finding the corre¬ 
sponding grain-t emperature distribution (GTD) f rom an observed 
SEP (see, e.g., T^ Ajoldsmith & Zhoul Il99ll : Ixie et al.] 1 19931 : 
lUi, Goldsmith & 7^ 1999h . Solving this inverse problem would 
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2 Mattsson et al. 


thus provide empirical constraints on the functional form of the 
GTD that can be used to improve dust estimates. In principle, there 
is a detailed correspondence between the dust SED and the GTD, 
but Ending the exact shape of the latter is unfortunat ely often a 
badly conditioned problem jHobson & Padmai]|l994h . Common 
practice is still to fit one or two single-temperature components to 
model the SED, which may affect the implied dust mass; in par¬ 
ticular, the mass of cold dust (~ 20 K) may be overestimated by 

pushing too much flux into this component. _ 

_ It has recently been suggested by, e.g., G12. [Richardson et al.l 

( l2013h and iTemim & Dwe^2013L henceforth TD13) that the treat¬ 
ment of grain temperatures may be important in the case of the Crab 
Nebula. Initially G12 used a two-temperature fit to the observed 
SED to derive dust masses of 0.25 Mq and 0.12M0 for silicates 
and carbon within the Crab respectively. In parallel with our work, 
TD13 have made a model of the EIR dust emission in the Crab 
Nebula explicitly taking the synchrotron radiation field and grain- 
sizes into account in calculating the radiation balance. They claim 
that the required dust mass to explain the SED is significantly re¬ 
duced in their model but this is due to a combination of the model 
as such and the use of a different set of optical constants (for which 
laboratory measurements are only available at T < 300 yum). We 
note, also, that deta iled modelling can give very different results 
compared to TD13 dOwen & Barlowl[2015ll . It is therefore worth 
elaborating on the effect of considering a cold-dust GTD and in 
particular, to consider how sensitive the dust masses obtained from 
SED fitting are to the assumed dust temperatures. We use a simple, 
computationally inexpensive, ‘model’ based on known constraints 
on the GTD, which can be applied in any environment regardless 
of whether the radiation field heating the dust can be specified or 
not. lKovacs et al.l (l2010h : lMagnelli et ^ d2012ll . have successfully 
modelled the GTD of star forming galaxies using a power-law in a 
similar manner. Attempts to solve the inverse problem mentioned 
above for molecular clouds have suggested t he GTD may be closer 
to an exponential form, however (see, e.g.. IXie et al.lll993h . The 
purpose of this paper is not to provide a precise and realistic model 
of dust emission, though. We merely aim to investigate how much 
of a difference introducing a range of grain temperatures makes in 
general, compared to other uncertainties in deriving the dust mass. 


2 FIR/SUB-MM EMISSION FROM DUST 

In this section we will just briefly summarise the physics underly¬ 
ing the conventional way of inferring dust masses and grain tem¬ 
peratures from EIR/sub-mm lluxefl 


2.1 Dust masses 

It is normally assumed that dust grains absorb and emit photons 
according to Kirchhoff’s law and that the source function can be 
described by a Planck function Bx\dAldv\ and that optically thin 
conditions apply in the EIR/sub-mm of the surrounding medium. 


* In the present paper we refer to ‘flux’ in the following way: the total flux 
F from a source is the energy output (luminosity) per unit surface area, e.g. 
in units of W m“^. The specific or monochromatic flux is per wavelength 
(or frequency) Fx is the same as ‘spectral flux density’ commonly used in 
observational work (usually with the unit being Jy = 10“^® W m“^ Hz“*). 
However, to keep the terminology simple we only use the term ‘flux’, and 
ask the reader to bear in mind that a subscript A means the flux at a specific 
wave length. 


If 2cm = 2abs. where 2 denotes the ratio of the effective and ge¬ 
ometric cross sections (Kirchhoff’s law) and pKxAx <K 1, where 
p is the gas density, kx is the opacity due to absorption and Ax is 
the geometric thickness of the surrounding medium (optically thin 
conditions), then the observed FIR/sub-mm flux is the sum of all 
reemitted flux from dust grains, where the flux from a grain of ra¬ 
dius a and temperature Td is 

ff(fl) = kD-^ Bx(Tx)a^Q,Ma)n(a)da. (1) 


In the equation above, D is the distance to the observer and n(a) 
is the grain-size distribution (GSD) by number normalised per unit 
volume. With the additional assumption that an ensemble of dust 
grains can be described with a single grain temperature Tj and that 
a/A <K 1 (the Rayleigh limit), the flux from this ensemble can be 
expressed as 

r-dust 3 Bd(rd) 2abs'‘8r 

" A - 

where 2ab.s('^) “ 2abs(d.a)/fl in the Rayleigh limit and V'gr is the 
total volume taken up by the dust. Multiplying both sides of the 
above equat ion with the bulk density of the grains p„, one ob tains 
the relation ( lHildebrandlll983l : rGalL Hiorth & Andersenll201 ll) 


3 2'b, MTi) " kxBxiTxY 


(3) 


where kx is a quantity to be referred to as ‘emissivity’ (or ‘absorp¬ 
tivity’, depending on the context) which has the same dimension as 
opacity, but reflects the optical properties of the grain material and 
should not be confused with the opacity of the surrounding medium 
Kx referred to above. 

Eq. lO conveniently evades the essentially unknown GSD. 
However, it is not obvious that the assumptions made to derive 
this simple relation between dust flux and mass are valid in gen¬ 
eral. In particular, it is only valid in the Rayleigh limit and if all 
dust grains have the same temperatures, where the latter implies 
that they should all have the same size as well. The Rayleigh limit 
only applies to emission at long wavelengths (absorption and scat¬ 
tering is mainly in the optical/UV), which may cause an implicit 
dependence on grain size with different grain temperatures as the 
result. In any realistic dust model one has to integrate over a suit¬ 
able size distribution and make use of a variety of dust compo¬ 
sitions/materials. That is, the dust mass corresponding to a given 
SED depends on the GTD W{Tx) = dMxjdTx. The flux from dust 
for such a distribution is given by 


W(Tx)Bx(Tx)dTx. (4) 

how 

We will in the following elaborate on how much of an effect W(Tx) 
has on the dust masses derived from FIR/-sub-mm SEDs. 


2.2 Grain temperatures 

The basic assumption underlying the relation between dust flux 
and dust mass is that it is the same dust grains that absorb light 
in the UV/optical or gain energy due to collisions with other par¬ 
ticles (usually electrons) that are also emitting radiation in the 
IR/FIR/sub-mm. The physical underpinnings for radiative and col- 
lisional heating are different, however, but the equilibrium temper¬ 
ature distributions for the grains may not be that different. 
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From flux to dust mass 3 


2.2.1 Radiative heating and cooling 


Radiative heating is normally due to absorption of radiation in the 
UV/optical and assuming the grains are in local thermal equilib¬ 
rium with the mean intensity of the radiation field we can equate 
the absorbed and emitted power, i.e. Pabs = Pem- The optical depth 
of the surrounding medium may in principle affect the energy ab¬ 
sorption but has little effect on the re-emission at long wavelengths. 
Thus, making use of (as previously defined) in the Rayleigh 
limit and assuming the surrounding medium is optically thin, we 
have 

r Q,^M,a)J^,,dA^AKa f Q:,^(A)B,(T,) dA, (5) 

Jo Jo 

where is the mean intensity of the radiation field. Note that 
2'^^ (the Rayleigh-limit approximation) can only be used on the 
right-hand side of Eq. which represents emission at long wave¬ 
lengths, and that Eq. © is strictly valid only locally (but can be 
generalised, resulting in an effective overall GTD). Conveniently, 
the emission at long wavelengths is approximately a power-law in 
A. Thus, replacing with QgCA/AoF^, integrating over wave¬ 
length and assuming that heating occurs in the grey-absorption 
limit (or large-particle limit) where 2abs is constant, we have then 
a simple power-law of the form 

/ a 

Pd(a) = Ps - . (6) 


where Tg and a., are scaling parameters. Note that even if one adopts 
an explicitly grain-size dependent 2abs, the grain temperature is still 
(locally) uniquely determined by the grain size under thermal equi¬ 
librium conditions. Eor further details on the above, see Appendix 

0 

The GTD VT(Pd) can be determined from 


W{Ti) = g>{a) 


da 


(7) 


where (fi( a) is the GSD in terms of mass. For a canonical MRN dis¬ 
tribution jMathis. Rumnl & Nordsiecklll977h (p(a) oc Thus, 
adopting an MRN distribution and the grey-absorption limit, the 
GTD is simply also a power-law W(Tii) oc Such a GTD 

has been recovered for hot dust (heat ed by short-waveleng th ra¬ 
diation) around active galactic nuclei ( IWang & Zhou|[l99^ . The 
grey-absorption limit may not be strictly applicable in other cases 
(e.g., SNRs) and the slope of the temperature distribution is there¬ 
fore likely steeper than in Eq. 0 and may deviate from the simple 
power-law form above, i.e., the scale temperature T, would in this 
case be a function of a. We note also that collisional heating should 
give rise to a GTD in a way similar to radiative heating (see Ap¬ 
pendix |B]l. 


2.2.2 Temperature fluctuations? 


Previous work has shown how the GTD is affected for small grains 
undergoing significant temperature fluctuations (non-equilibrium 
conditions) as a result of being h i t by energetic photons or 


other particles (see. e.g.. PurcelJ 19761; 

Aannestad & Kenvonll 19791; 

Draine & Andersonll 19851; iDwekll 1986 

; [Draine & Li[[200l[; [Draind 

2003|. and references therein). The heat capacitv determines how 


much energy the dust particle can hold, i.e., the energy in an ab¬ 
sorbed photon will partly heat up the material and the remainder 
will be emitted at longer wavelengths. The amount of energy a grain 
can hold is proportional to its mass, which is why small grains do 


not easily obtain an equilibrium temperature if heated with short- 
wavelength radiation. For radiative heating, we have that an en¬ 
ergy increment dE is causing a temperature increment according to 
dE = Cv dTi, where Cv is the heat capacity at a constant volume. 
The heat capacity for a spherical grain is, in the low-temperature 
limit of Debye’s model (TgITo <sc 1, with Tjg the Debye tempera¬ 
ture). 



16;r^ kpg,a^ I Tg ? 
5 rax \Td/ 


( 8 ) 


where N is the number of monomers in the grain, mx their mass 
lAannestad & Kenv^l 19791 : IPraine & Andersoiilll985b . The tem¬ 
perature increase is effectively ATg ~ 2abs(‘*> T) Cy' iiphot ~ 
where iiphot is the energy of the photon and gabs is the ratio of the 
effective and geometrical cross sections of the grain. Inserting suit¬ 
able numbers shows that a single UV photon could suffice to raise 
the temperature of a cold (~ 30 K) nano-sized grain almost an order 
of magnitude. Very small grains will usually also cool rapidly and 
regain their i nitial temperatures of ~ 30 K in t ypically a couple of 
hours or less l lPurce Iilll976l:lDraine & Lill200lh . Thus, at any given 
time, a population of cold, very small grains may therefore have 
wide range of temperatures, leading to a GTD even for a popula¬ 
tion of grains with exactly the same sizes. 

At first glance it may seem like the above factors compli¬ 
cate the picture a lot. But, fortunately, the effect on the estimated 
mass of cold dust is negligible. First, most of the mass is ex¬ 
pected to be in the large grains (cf. the MRN size distribution 
in terms of dust mass, dpi/da = ip(a) oc and dust in 

SNe appears to have a stron g bias towards large grains in general 
according to recent re sults jGall et alj|2014l : IWesson et alj|2014| ; 
lOwen & Barlowll2015ll . Not more than a few percent of the mass 
can be in gr ains small enough to undergo significant t emperature 
fluctuations jPurcelllfl9761 : lAannestad & Kenvon|[l979h . From eq. 
0 we can also see that size can easily compensate for the fact 
that the grains are cold. Second, if temperature fluctuates are im- 
portan t, they will cause e mission in the infrared and near infrared 
bands l lLi & Drainell200lll . Cold dust radiates almost exclusively in 
the FIR/sbmm bands. Thus, in conclusion, we do not have to worry 
greatly about temperature fluctuations as long as we are dealing 
with cold dust and a steep GTD. 


3 THE EFFECTS A GRAIN-TEMPERATURE 
DISTRIBUTION 

Under most circumstances it is inevitable that there is a distribution 
of dust temperatures rather than distinct representative dust tem¬ 
peratures for specific dust components (e.g., a cold, warm and hot 
component) as we have touched upon above. This fact has pre- 
viously been pointed out as a caveat (see, e.g., G12, TD13 and 
[Richardson et alj2013l . in the case of the Crab Nebula). Such a dis¬ 
tribution must be very steep and narrow not to create SEDs which 
are inconsistent with the featureless ‘bumps’ associated with cold 
dust, but it can still have a significant effect on the derived dust 
mass. In this Section we describe a simple multi-temperature model 
for the FIR/sub-mm dust emission. We use the Crab Nebula as a 
test case, to show that an observed SED can be reproduced with 
very different dust masses, and then continue with a more general 
analysis of how multi-temperature SED fits may differ from two- 
temperature fits, by generating a large mock sample of SEDs from 
simple GTDs. 
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4 Mattsson et al. 


3.1 Grain-size or temperature distribution? 


There are essentially just two ways to incorporate the effect of a 
range of grain temperatures when modelling an SED: either one 
has to (computationally) find the GTD from an assumed GSD and 
a (possibly) known radiation field, or, one could make a direct as¬ 
sumption or estimate of the GTD (or its functional form, more pre¬ 
cisely). TD13 preferred the first approach, while we will here ex¬ 
plore the latter. 

There is no actual advantage of any of these two approaches 
over the other in terms of their physical correctness (we simply 
have too little information to distinguish between the two), though 
it is clear that an assumption directly regarding the GTD is much 
simpler to deal with. There is also a direct link between the ob¬ 
served SED and the corresponding GT D that could (in theory) 
be used to constrain the functional form l Xie. Goldsmith & Zhoul 


I 1991 I : IXie et al.lll993l : iLi. Goldsmith & Xid 19991) . In principle ane 

could thus construct a fairly consistent model even without any 
knowledge about the heating radiation field. 

One could of course argue that the best option would be to 
use a GSD that results from theoretical modelling of SN dust for¬ 
mation and then chose an approach similar to that used by TD13. 
But unfortunately th e models do not offer a consistent picture. 
iNozawa et al.l (l2003h suggest the GSD produced by SNe may be 
somewhat flatter [n(a) oc except for very large grains, where 

the MRN distribution seems to be recovered. But the effective GSD 
mayb e even flatter still, if small grains are destroyed dNozawa et al.l 
l2007h . and the upper and lower size limits do not correspond to 
sharp cut-offs, so the assumption of a power-law GSD made by 
TD13 is not obviously the optimal ansatz. 


3.2 Multi-temperature SED 

Modelling of the dust contribution to the SED is often done by fit¬ 
ting one or two components (cold and warm dust), with weighting 
factors that specify the mass contribution from each component. 
We generalise this to N components, but with the weighting factors 
constrained by a GTD assumed to follow a power law or expo¬ 
nential as motivated above. Thus, we maintain the same number of 
parameters (four) when fitting the SED, while at the same time dust 
grains can have an arbitrary range of temperatures. 

Regarding the functional form of the GTD WiT^), the energy 
balance (Eq. |5j suggest that if a grain has a temperature Ti and 
another grain has a temperature T 2 > Ti, then grains with tem¬ 
perature T 2 cannot be larger than grains with temperature Ti (see 
also the Appendix |3. That is, dWIdT^ < 0 for all T^. Eurther- 
more, because 2abs, at any given wavelength, normally is a smooth 
function of a without very prominent features, there is no reason 
to expect that W{Ti) is a complicated function of T^. It is there¬ 
fore reasonable to assume that the GTD is a smooth, monotonously 
decreasing function of grain temperature. Ste ep power-laws have 
previously been considered in other con texts dAguirre et al.ll2003l : 
ikovacs et al.ll201(]l:lMagnelli et alj201 j) . 


3.2.1 Power-law 

As discussed previously, radiative heating that occurs in the grey- 
absorption limit corresponds to a power-law GTD W{Ti) where the 
GSD is a power-law. If the dust is heated by short-wavelength radi¬ 
ation (UV, X-ray etc.) this is in fact a good approximation of reality. 
We have therefore as good reasons to try a power-law form as our 


ansatz for the effective GTD as there are reasons to expect a power- 
law GSD. 

With W{Ti) - Wo (TilToy”, where we set Tq = we can 
model the flux from an ensemble of dust grains with different sizes 
and temperature using a function 


S.fThigh, Tiow; a) = So(N) ^ B,(r,) 


(9) 


in which T, = i A^'fEhigh - Tiow) + Tiow, where Thigh and Tiow are 
the high and low temperature cut-off, respectively. So is a constant 
such that S a represents the blackbody flux from all N components. 
This function replaces the Planck functions and the weighting co¬ 
efficients in, e.g., a two-component model and the exponent a can 
be treated as a free parameter. If N is large (we assume N = 100000 
in subsequent applications), the fitting result corresponds to a con¬ 
tinuous GTD W{Ta) as in Eq. dUl. 


3.2.2 Exponential 

The power-law form of W(rd) is strictly valid only under very spe¬ 
cial circumstances, as we have pointed out in Section 12.2.1 1 It is 
more realistic that the slope of the G TD changes with the grain 
temperature Tj. The GTDs derived by IXie et^ ( Il993h from the 
SEDs of cold molecular clouds are very close to exponential for 
Td <: 25 K and detailed computer simulations of molecular clouds 
do in fact suggest this may be a realistic functiona l form of the 
GTD of a single dust species llBethell et al.l |2003) . If an expo¬ 
nential form of IV(rd) works well for molecular clouds it seems 
reasonable it could work as well for cold dust in SNRs. Thus, 
we try also an exponential model. Using an exponential GTD, 
W(Ti) = W^o exp(-7'd/To), we have in analogy with the power law 
case, 

JV-l 

SdIThigh, Plow; To) = 5o(Al) 'Yj 

i=0 

where To is a free scaling parameter. Since a GTD must be steep, a 
realistic value of To is only a few K. 

Power-law models used on very c old dust suggest the GTD is 
flattening out at low dust temperatures jAguirre et al.ll200^ . which 
further motivates that we should consider a GTD of the exponen¬ 
tial form above (see Pig. [T] for an example of how this flattening 
naturally occurs at low Td for an exponential GTD). 


Bd(T,), 


( 10 ) 


3.2.3 Temperature cut-offs 

Obviously, there has to be a lowest and a highest possible grain tem¬ 
perature, because grains are neither arbitrarily small, nor arbitrarily 
large (see AppendixlAt. We use these cuts on the temperature distri¬ 
bution as fitting parameters as well (see Eq.|^, which corresponds 
to the integration limits above. Initial values for Ti^,, and Thigh may 
be determined by a two-temperature fit, which stabilises the fit¬ 
ting procedure. One may argue that the lower limit Tio* may also 
be fixed to this value, since the lowest possible dust temperature 
is essentially defined by the shape of the PIR/-sub-mm tail of the 
SED (see fig. 4 in G12). But this is not necessarily a good approach, 
since the flux contribution from very cold (and large) grains may be 
small, while their contribution to the dust mass is significant. Note, 
also, that this is related to the upper grain-size limit as well (the 
coldest grains are the largest ones), while the upper temperature 
limit, and thus also the lower grain-size limit, is harder to constrain 
(we will return to this later). TD13 assume an upper grain-size limit 
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Figure 1. Grain temperature distributions for the dust SED models with 
Flow fixed to the temperature of the cold component in a two-temperature 
fit of the Crab Nebula. The temperature interval used to model the SED for 
each model is shown by solid lines. See also Table[T] 


which, according to their model, corresponds to approximately the 
temperature of the cold component in a two-temperature fit. There¬ 
fore, in our ‘test application’ to the Crab Nebula below, we have 
included a case where Tiow is fixed to the temperature of the cold 
component in the two-temperature fit, in order to compare with the 
results of TD13. 


3.3 Test case: fitting the Crab Nebula SED 


We have chosen to use the Crab Nebula to test our simple multi¬ 
temperature component model since it has a wide FIR/sub-mm 
SED which seems to suggest a range of grain temperatures is possi¬ 
ble (see also the papers by TD13 and lRichardson et al.r2013l where 
a range of grain temperatures is discussed as well). SN 1987A, on 
the other hand, has a very narrow SED which appears consistent 
with a singl e-temperature populat ion with an extremely low dust 
temperature dMatsuura et abll^l ll) and thus makes it irrelevant as 
a test case in the present context. 

There is also a third well-observed remnant that has received 
a lot attention beca use of its seemingly large dust mass, Cas A (see 
IPunne et al.ll2009l and references therein), but Cas A is known to 
be a complicated case. It is far from established that the FIR/sub- 
mm part of the SED primarily reflects dust in the SNR, because 
of foreground contamination from a sp i ral arm i Dunne et'^l2003l : 
iKrause et alj|2004l : iGomez et al.ll200^ : ISarlow et alj|20l3r Thus, 
we decided to not consider Cas A as a test case either, although the 
shape of the SED would suggest a fairly wide temperature distribu¬ 
tion as seen in the Crab Nebula. 

We use the photometric data from G12 to plot the SED 
of the Crab Nebula (Eigure with wavelengths ranging from 
1 -1000 yum. The Crab Nebula is a somewhat special object in that it 
has a strong synchrotron radiation field originating from the pulsar 
wind nebula (PWN). The PWN is the main heating source, a fact 
that TD13 take advantage of in their model (i.e., in this particular 
case, the heating source can also be specified quantitatively, but this 
is not the case for other SNRs). Eor our purposes, we only have to 
subtract the synchrotron component from the SED, which is rela¬ 
tively easy since it is well-described by a power-law (see G12 and 
references therein). The integrated fluxes need also to be corrected 


for line emission. In some cases it is small, e.g., only 8.7% of the 
total flux at 100yum and 4.9% of the 70yum (see Table 2 in G12). 
But the 24 yum flux due to dust emi ssion is 43% of the synchrotron 
subtracted flux jTemim et al.l201 j) . which is important to take into 
account. 

Since G12 obtained a valid fit wit h a canonical two- 
temperature fit with astrono mical silicates dPraine &J^ej jjj84|: 
Weingartner & Drainell200i]) and amorphous carbon l IZubko et al.l 


1996h . we first modelled the SED using astronomical silicates but 
now with the two forms of GTDs described in Section[3 The range 
of dust temperatures (Table [T) obtained through the SED fitting 
corresponds to a steep distribution favouring cold dust (see Eig. 
[TJ. Since the low-temperature limit Tiow is decisive for the in¬ 
ferred mass, we have considered two cases: (1) a lower limit on 
the dust temperature set to the cold dust component temperature 
derived in G12 (again, see Table [T) and (2) a lower limit which 
is treated as free parameter, albeit with an absolute lower limit set 
by the cosmic background temperature Tcmbr = 2.73 K. We will 
later discuss the effects of assumptions about Tiow in more detail 
(Section 3.4). The upper limit is treated as a free parameter for all 
cases, but is unsurprisingly rather similar to the temperature of the 
warm component as obtained from a two-temperature fit (see Ta¬ 
ble [T]l. The models with a range of temperatures (full-drawn black 
and dot-dashed red curves in Eigs.|2and[3 are as good fits to the 
data (lower;y^) as the two-temperature fits (the blue dashed lines). 
But the corresponding dust masses are quite different. Assuming 
silicate dust and Tiow fixed at the value obtained from the two tem¬ 
perature fits, only ~50% of the dust mass is required (compared 
to the two -temperature fit). U sing amorphous carbon grains (data 
taken from IZubko et al.ll996h with a range of temperatures instead 
of silicates yields similar results, though with a revised dust mass of 
60-70% compared to the carbon grain model in G12. The second 
case, where Tiow is a fitting parameter, leads to lower Tiow values 
and thus significantly higher dust masses. With a power-law GTD 
the dust masses, assuming silicates as well as amorphous carbon, 
are roughly doubled compared to the two-temperature fits. Using 
an exponential GTD, the fitting algorithm pushes Tiow to the min¬ 
imum value Flow = Fcmbr = 2.73 K. Maximising the amount of 
very cold grains like this suggests dust masses which are a factor 4- 
6 higher than those obtained in the two-temperature fits. Obviously, 
this last result is not very realistic, but it clearly demonstrates why 
fitting simple SED models to data can be dangerous. 

The power-index values we obtain for t he power-law model s 
agree well with the v alues obtained by, e.g., iKovacs et al.l (l2010l) : 
iMagnelli et al.l (1201 2h who used a similar GTD app roach for de¬ 
riving the dust mass in galaxies. iKovacs et al.l (1201 OL see also ref¬ 
erences therein) discussed that a power index in the range a = 
6.5...8.5 is expected in diffuse media, while for dense interstellar 
media a = 5...1 is more likely. The upper end of these ranges cor¬ 
respond to an effective emissivity index p = 2, which is appro¬ 
priate for, e.g., astronomical silicates. Thus, a SNR (which can be 
regarded as a dense medium) with silicate dust should have ff = 7, 
which is exactly the value we have obtained for the Crab Nebula 
with Flow fixed (see Table [TJ. For amorphous carbon dust, which 
has p x: 1, we should expect a k 6 , in agreement with our results 
in Tabled] 

In principle, the dust masses obtained with a fixed Flow could 
lower the tension between the expected metal budget and the 
amount of dust formed in the ejecta. The mass of the progenitor 
star is constrained by the fact that one can put relatively strict con¬ 
straints on the mass of the neutron star, M„s :=! IAMq, which sug¬ 
gest a progenitor mass below 13 Mq according to theoretical mod- 
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Table 1. Comparison of parameters derived from fitting the SED of the Crab Nebula with two-temperature and GTD models. ^low and 7’high denote the low 
and high cut offs of the GTD as well as the temperatures of the cold and warm components in the two-component models, a is the resultant power-law index 
for the GTD. 


Model 

Dust type 

Pgr 

[g/cm^] 

Flow 

IK] 

Thigh 

[K] 

To 

[K] 

a 

Mi 

[Mq\ 


Remark 

A 

Astron. silicate 

3.3 

28.1 

55.6 

- 

. 

0.25 

8.54 

Two-temperature fit. 

B 

Astron. silicate 

3.3 

28.1 

72.6 

- 

7.0 

0.14 

7.98 

Power-law GTD, Aow fixed. 

C 

Astron. silicate 

3.3 

17.4 

68.7 

- 

5.2 

0.51 

7.98 

Power-law GTD, Tiow bee. 

D 

Astron. silicate 

3.3 

28.1 

78.3 

6.38 

- 

0.12 

8.23 

Exponential GTD, fixed. 

E 

Astron. silicate 

3.3 

2.73 

74.1 

6.37 

- 

1.14 

7.87 

Exponential GTD, Tiow bee. 

F 

Amorphous carbon 

1.81 

33.8 

63.4 

- 

- 

0.12 

8.19 

Two-temperature fit. 

G 

Amorphous carbon 

1.81 

33.8 

78.7 

- 

6.5 

0.077 

7.66 

Power-law GTD, Tiow fixed. 

H 

Amorphous caibon 

1.81 

19.8 

74.0 

- 

4.1 

0.22 

7.79 

Power-law GTD, Tiow bee. 

I 

Amorphous caibon 

1.81 

33.8 

83.9 

7.79 

- 

0.072 

7.75 

Exponential GTD, Tiow fixed. 

J 

Amorphous carbon 

1.81 

2.73 

76.2 

11.0 

- 

0.71 

7.73 

Exponential GTD, bee. 






Figure 2. SED model fits to the Crab Nebula using a power-law GTD 
with photometric data from G12. The dashed blue lines correspond to the 
two-component models by G12, while the solid black and read dot-dashed 
lines show temperature-distribution models with Tiow fixed and as a firing 
parameter, respectively. The thin dashed black line shows the estimated syn¬ 
chrotron component. The upper panel show models based on silicate dust, 
while the lower panel shows models based on amorphous carbon dust. 


els (see, e.g.. lWooslev & Weaveilfl99^ . But we note also that the 
chemical abundances seem to be in better agreemen t with a pro¬ 
genitor above 11 Mq jMacAlpine & Satterfiel j|2Q^ . In general, 
it seems a progenitor of relatively low mass i s favoured, which 
is supported by the slow expansion velocity i Fesen et al. ira; 
ISmithiHool) . A recent assessment by ISmithI (12013 ) has strongly 
constrained the progenitor mass to 8-10 Mq, suggesting a super- 
AGB star that suffered an electron-capture SN rather than an Fe 
core collapse event, which puts a severe limit on the metal bud¬ 
get. The masses found by G12 are then only marginally consis¬ 
tent w ith the metal budget i f one accepts the nucleosynthetic mod¬ 
els by IWooslev & Weav^ d 19951) . A 13 Mq star in such case al¬ 


Figure 3. Same as Fig.|2]but for an exponential GTD. 


lows 0.37 Mq of silicates (assuming an effective mass number 
^sii = 170) and 0.11 Mq carbon dust to be formed. Thus, while the 
dust masses found by G12 are not obviously overestimates, a model 
suggesting close to 100% dust-condensation efficiency is not very 
convincing since one expects a non-negligible sublimation rate in a 
SNR. The lower dust masses may seem like a reasonable and con¬ 
servative choice, but one has to be aware that the total SED cannot 
provide conclusive evidence. 


3.4 How much difference does a multi-temperature fit make? 

3.4.1 Shape of the SED 

The range of dust temperatures is directly connected to the width 
of FIR/sub-mm (dust) bump in the SED. ‘Warmer’ grains (with 
temperatures above the coldest grains, that is) will inevitably add 
flux on the short-wavelength side of the dust SED. What this extra 
flux will look like in the SED depends on the dust type, though. 
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Table 2. Pai ameter ranges for the grid of artificial SEDs. 


Parameter 

Range 

Step size 

Flow 

20-50K 

0.3 K 

Fhigh 

60 - lOOK 

0.6 K 

To 

3-9K 

0.06 K 

a 

4-9 

0.05 


Astronomical silicates can provide a tell-tale signature - the lOyum 
feature - that reveals the presence of warm dust (T^ >100 K), while 
the carhon-dust FIR/sub-mm SED is mostly featureless. The 10 /jm 
feature cannot arise from cold dust, however. But the location of 
the peak wavelength and slope of the short-wavelength tail poses a 
constraint on the cold-dust GTD regardless of the presence of any 
dust-emission features. 

The kurtosis (‘peakedness’) of the dust SED can also be af¬ 
fected by the GTD. In particular, using a continuous GTD may re¬ 
sult in more flux in the middle of the SED than in the case of a 
two-temperature model. This is clearly seen in the fits to the Crab 
Nebula with a fixed low-temperature limit (see Section Figs.j^ 
and[3. Thus, a continuous GTD may differ slightly from the two- 
temperature model in terms of the shape of the resultant SED. This 
affects the required dust mass of the model because the fit to the 
observed SED will be different. An increased kurtosis in the model 
SED should result in a lower dust mass, which is also what we ob¬ 
tained in our application to the Crab Nebula. 


3.4.2 Reducing the dust mass? 


In the continuous case (the number of temperature components 
N oo), we can write 5 a as a temperature mean of the Planck 
function Bx weighted by the GTD W{Ti), 


Sx 



W(Tx)dTx 


rAigh 


W(Ti)Bx(Ti)dTi, 


( 11 ) 


which replaces Bx in the derivation of the dust mass from a given 
SED. The total dust mass obtained from a two-temperature fit is 
completely dominated by the cold component (see, e.g., G12) and, 
as we have pointed out above, it seems reasonable to associate the 
low temperature cut-off Tiow with the temperature of the cold com¬ 
ponent Tcoid from the two-temperature fit (although this is not ob¬ 
viously the case, which we will discuss later). Hence, if we com¬ 
pare Sx(Tiov, Thigh', To) and BxiTio„), we would have an approxi¬ 
mate measure of how much the GTD affects the implied dust mass 
for a given set of parameters, if the cold-dust temperature is equal to 
the lower temperature limit Tiow of the GTD. Since we are consid¬ 
ering thermal radiation at long wavelengths, we may also approxi¬ 
mate the Planck function using the fact that e*^'-l ~ xforx» l.It 
is then straight forward to show that, for W{Ti) = Bo expf-Td/ro), 
Tq <k Thigh (very steep GTD) we have at large wavelengths 


Bx Flow 


( 12 ) 


which is because the effect of a wide range of dust temperatures 
is in such case limited and depends on the balance between cold 
and warm dust. That is, the effect of a realistic GTD (which must 
be steep) on the implied dust mass depends mostly on the temper¬ 
ature of the coldest dust. Moreover, using Eq. J12t and since Ti^^ 
is typically a few times Fq (if the temperature Ti^^ is that of the 
two-temperature fit), one can easily verify that the increase in flux 


From flux to dust mass 


from a GTD compared to the two-temperature fit is not likely more 
than 20 - 30%. 

We have also computed Eq. (ID by numerical integration us¬ 
ing an exponential as w ell as an power-law GTD w ith Qabs for 
‘astronomical silicates’ (Weing ptner & Drainel 1200 ih and amor¬ 
phous carbon (Zubko et all 19961) . Because this integration is com¬ 
putationally fairly inexpensive, we generated a large, dense grid 
of artificial SEDs with either a fixed lower temperature limit or 
a fixed upper temperature limit and a range of values for the re¬ 
maining parameters (see Table |D- For each SED we then fitted 
two-temperature models with the temperature of the cold compo¬ 
nent fixed to Fcoid = Flow and with Fcoid as a free parameter. We 
then compared the required dust masses, i.e., we computed the ra¬ 
tio These ratios are shown in Figs.|4l[3|6]and|7] In 

total we computed 4 • 10^ SEDs and made two-temperature fits all 
of these, but we only show a subset of 80000 SEDs here (since the 
results for silicates are very similar to those for amorphous carbon 
and we do not display all the various cases of fixed upper and lower 
temperature limits either). 

As can be seen on Figs.|4]and|5l the -ratios sug¬ 

gests the introduction of a GTD (with Fi^w = Tcoid) only lowers 
the implied dust mass by ~ 30 - 40% for realistic values of Flow, 
Thigh and To in case of an exponential GTD (left panels in Fig. |4] 
and|D and a similar result is obtained for the power-law case (right 
panels in Fig.|4]and[D- Comparing the Figures, it is also clear that 
the lower temperature limit Flow is more important than the upper 
limit Flow for any reasonable GTD (which must have To Fhigh or 
a ~ 5...8), as we predicted from Eq. I ll2l l above. 

The ~ 50% (or 30 - 40%) reduction of the (silicate or carbon) 
dust mass that we obtained by including a GTD in the SED mod¬ 
elling of the Crab Nebula is in fact due to two factors that make 
comparable contributions: the increased flux due to the addition of 
a range of grain temperatures (higher average Fd) and a mere dif¬ 
ference in the fit due to the slightly different shape of the model 
SED. The latter can be seen in Figs. |2] and as well as Fig. 5 in 
TD13. The combined effect amounts to & factor of two for silicates 
in the Crab Nebula, which is likely the largest reduction one can ex¬ 
pect in employing a GTD in general (not only for the Crab Nebula). 
This is expected since in case of an exact fit to the SED, the rela¬ 
tive mass reduction would be the same as the relative flux increase 
from adding more warmer dust (typically 20 - 30%). As is evident 
from figure|4]and|5] the expected dust-mass ratio, comparing a two- 
temperature fit with a GTD-model fit, would only in extreme cases 
(e.g., very cold dust) reach above two. 

It is worth stressing that the mass reduction we describe above 
is totally dependent on the assumption that the temperature of the 
cold-dust component in a two-temperature fit is the same as the 
lower temperature limit of the GTD (Flow = Fcoid)- This assumption 
is similar to assuming an upper-limit grain size, as in TD13, and has 
essentially the same effect: it limits the amount of dust mass. How¬ 
ever, we performed further numerical experiments which showed 
that if one allows the GTD to reach down to very low dust temper¬ 
atures, a two-temperature fit to that GTD may actually predict the 
same or even a higher dust mass, compared to the two-temperature 
fit by G12. The effect of the extra flux from warmer dust on the pre¬ 
dicted dust mass may, in reality, be counteracted by the presence of 
dust colder than the cold component of a two-temperature fit. Such 
cold grains may affect the dust mass without making a significant 
contribution to the SED, thus still yielding a good model fit. 

In Figs.|^and|7]we show the implied dust masses from a two- 
temperature fit to a grid of SEDs generated from simple GTDs com¬ 
pared to the ‘true’ dust masses corresponding to the adopted GTDs 
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Figure 4. The implied dust masses from a two-temperature fit to a grid of SEDs generated from simple GTDs compared to the ‘true’ dust masses corresponding 
to the adopted GTDs. Left panel: the effect on the dust mass assuming an exponential GTD with Tygi, = 100 K and To the scaling temperature of the GTD. 
Right panel: same as left panel, but for a power-law GTD with power index a. The temperature of the cold component in the two-temperature fits are assumed 
to be the same as the low-temperature cut-off in the GTD, which represents a reasonable upper limit to the grain size. These figures may not display well in 
b/w, but are available in colour in the online version of the journal. 


as in Figs. |4] and 12 except that the temperature of the cold com¬ 
ponent in the two-temperature fits is here treated as a free fitting 
parametej^ The dust masses inferred from the two-temperature 
fits are typically ~ 10 - 40% lower than the ‘true’ dust masses. 
The cold component temperature is typically 5 - 10 K higher than 
the low-temperature cut-off in the GTD, although there is no gen¬ 
eral scaling relation (see the contours in Figs. |2 and [7]l. With a 
fixed high-temperature cut-off at Thigh = 100 K the low-temperature 
range Tiow = 25 - 35 K represents a special case: for an exponential 
GTD the two-temperature dust mass has a local minimum, and for 
a power-law GTD there is local maximum in this interval. This 
special low-temperature range depends somewhat on the choice of 
Thigh and for narrow GTDs such a temperature range may not exist. 
But it coincides with the range of dust temperatures that is the most 
interesting for cold dust in supernovae. It seems, therefore, that the 
effect on the inferred dust mass depends on the form of the GTD, 
which indicate that detailed radiative transfer models may he the 
best option as it will provide the best possible information about 
the GTD. Overall, our numerical experiments presented here, sug¬ 
gest that a two-temperature fit may not be as bad as it may seem at 
first, judging from the results of TD13 as well as our results pre¬ 
sented in Figs.[2]and[2 but simple SED fits are still not reliable. 


^ Note that this is different from the Crab Nebula fits with as a free 
parameter, although it illustrates the same phenomenon. 


3.5 Temperature-size relation 

In an equilibrium model, there is a direct relation between grain 
size and grain temperature for a given (invariant) radiation field, 
which should be evident from Section [2.2.1l (but see also Appendix 
\3- In a simple model, like the one we have used here, it is not 
meaningful to discuss grain sizes in quantitative terms based on the 
steepness of the GTD and temperature cuts as obtained from SED 
fitting. The reason is that Eq. 0 is a local relation and adopted 
GTD is a global approximation. But locally the grain size must 
uniquely determine the grain temperature as long as nothing breaks 
the equilibrium. The steepness of the GTD is also directly depen¬ 
dent of dTi/da, i.e., the temperature-size relation is fundamental to 
the connection between the dust mass and the dust SED (see Eq. 
[7] and Appendix IATi. But one should also remember that there is a 
degeneracy between the upper size limit of the grains and the slope 
of the GSD, which can beco me problematic even in models with 
detailed radiative transfer (see lOwen & Barl^l2015h . 

The upper temperature limit (corresponding to lower size 
limit) may have bearing on the implied dust mass obtained from the 
SED (see Eigs.|2and|3 since the small grains usually only make up 
a small fraction of the dust mass, but still contribute significantly 
to the SED since they are slightly warmer than the large grains (the 
flux from a grain is essentially proportional to the temperature). 
TD13, on the other hand, argue that the choice of Amin has no par¬ 
ticular effect on their results for the Crab Nebula (according to their 
Monte Carlo simulation) which is a result we can only reproduce 
when Flow is fixed and the GTD has a certain slope. Possibly, this 
is related to the fact that their model yields dT^jda k: Q for small 
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Figure 5. Same as Fig.|4]but with the lower temperature limit fixed at Thigh = 25 K and variation of the upper limit. Left panel shows the case of an exponential 
GTD and the right panel a power-law. 



Figure 6. Same as Fig. [4| but without assuming that the temperature of the cold component in the two-temperature fits is the same as the low-temperature 
cut-off in the GTD. 
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Th,Bh (K) (K) 



0.68 0.70 0.72 0.74 0.76 0.78 0.80 0.82 


0.74 0.76 0.78 0.80 0.82 0.84 0.86 


Figure 7. Same as Fig. |4] but without assuming that the temperature of the 
cut-off in the GTD. 

grain radii a, though it is unclear why becomes size independent 
for small grains in their model, which is also an equilibrium model. 


4 DISCUSSION AND CONCLUSION 

Since there is known to be a distribution of dust-grain sizes, there 
must also be a distribution of grain temperatures (GTD) - even in 
case of thermal equilibrium, which should apply to cold dust. We 
illustrated how this can increase the emission from any type of dust 
component, whether it is silicates, carbonaceous dust or another 
composition. We applied a GTD model to the SED of the Crab 
Nebula, which can be explained using thermal emission from solid 
dust grains at a range of dust temperatures (rather than the canon¬ 
ical two-temperature component model). A range of dust temper¬ 
atures lower the required dust mass by ~50% and 30-40% for as¬ 
tronomical silicates and amorphous carbon grains compared to re¬ 
cently published values (0.25Mq —> 0. 14Mq and Q.Q\2Mq —> 
0.0085Mq, respectively), but the implied dust mass can also in¬ 
crease by as much as almost a factor of six (0.25Mq —» 1.14Mq 
and 0.12Mq —> 0.7IMq) depending on assumptions regarding the 
sizes/temperatures of the coldest grains. In general, we find/confirm 
that: 

(i) The width of the SED determines how much small warm 
grains can contribute to the EIR/sub-mm flux excess. Only SNRs 
with relatively wide SEDs may show GTD effects worth consider¬ 
ing. 

(ii) Introducing a GTD has a significant but limited effect on 
the derived dust mass compared to a two-temperature fit. With the 
lower temperature limit of the GTD set equal to temperature of the 
cold component of the two-temperature fit, the implied dust mass 


component in the two-temperature fits is the same as the low-temperature 


is typically ~ 50% larger for the two-temperature fit. But without 
this coupling there is not necessarily any dust-mass reduction 
resulting from the introduction of a GTD (the inferred dust mass 
may in fact increase). 

(iii) The difference in shape of the SED (e.g., its ‘peakedness’) 
for a two-temperature model compered to a GTD model, in com¬ 
bination with the flux uncertainties, lead to different fitting results, 
which may affect the implied dust mass as much as the extra flux 
from warmer grains added due to the GTD. 

We have therefore shown that introducing a GTD may predict 
a dilferent dust mass than a canonical two-component model. But 
we also demonstrate that the effect is limited: we have shown that 
this depends strongly on the temperature of the coldest dust and 
how well-constrained the observed SED is in the EIR/sub-mm. It is 
not appropriate to claim that introducing a GTD (or a grain-heating 
model as in TD13) gives us better constraints on the dust mass of a 
SNR, such as the Crab Nebula. Eirst, we cannot know whether there 
should be a low-temperature limit in the GTD with a value similar 
to the temperature of the cold component in a two-temperature fit, 
even if this seems a reasonable assumption. The dust mass is uncer¬ 
tain by at most a factor of a few due to this. Second, the distances 
to many Galactic SNRs are relatively uncertain. The Crab Nebula 
is no exception: G12 adopt D = 2 kpc, which is exactly in the mid - 
dle of the range D = 1.5 - 2.5 kpc given bv lKaplan et alj (l2008h . 
Assuming D = 2.0 ± 0.5 kpc the uncertainty of the dust mass is al¬ 
most a factor of three. Third, uncertainties in the optical and struc¬ 
tural properties of the dust component amounts to at least a factor 
of a few: the emissivity of dust may vary considerably according 
to some observatio nal estimates (see, e.g.. I Alton et alfcOOOL 12004 : 
IPasvra et al. l2005h and if volatiles (e.g. ice mantles, which we will 


cold 
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discuss in a forthcoming paper) makes up a significant part of the 
dust, there may he an additional factor of two in the uncertainty. 
Combining all of the above, the total (maximum) uncertainty range 
spans at least an order of magnitude, somewhat depending on avail¬ 
able constraints. A major point is that while the present GTD ap¬ 
proach is not more accurate than the two-temperature approach, it 
demonstrates that simple SED fitting cannot really constrain the 
dust mass. 

One may then ask whether the GTD is ultimately an improve¬ 
ment. Here we would argue that incorporating a more physical 
model for dust heating and FIR/sub-mm emission from dust grains 
is indeed always a qualitative improvement, but the most concern¬ 
ing uncertainty in the model - what the dust is actually made of - 
still remains. However, we have demonstrated that the GTD does 
not always reduce the dust mass derived from SEDs (as implied by 
TD13), but also that the GTD-related uncertainty is typically not a 
dominant source of uncertainty i.e. the effect of grain temperatures 
is usually a part of the total uncertainty in deriving dust masses 
from the SED. 

The degree of dust condensation (the fraction of condensible 
material that end up in dust grains) is unlikely close to 100%, but 
the uncertainties in the conversion from FIR/sub-mm flux to dust 
mass and the amount of metals available for dust formation prevent 
precise estimates. In conclusion, the efficiency of dust production 
in supemovae remains poorly constrained, even if we would con¬ 
struct a sophisticated model of dust emission. This does of course 
not mean that radiative transfer models based on dust-grain pop¬ 
ulations with a range of grain sizes and temperatures are not im¬ 
portant (and needed) tools for converting infrared to submm SEDs 
to dust masses. We conclude that one should avoid using simple 
SED fits as far as possible whenever there is evidence of a range 
of grain temperatures. But the overall error in the amount of dust 
formed in supernova is still largely due to the uncertainties in dust 
composition, structure and behaviour of optical constants and thus 
determining the mass of dust remains extremely difficult. 
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APPENDIX A: HEATING, COOLING AND EQUILIBRIUM 


If the dominant source of heating is radiation (typically in the 
UV/optical), the dust grains may be regarded as being in local 
thermal equilibrium with the mean intensity of the radiation field 
around them, i.e., for grains of a specific radius a and temperature 
Pd we have an energy gain (heating) due to absorption given by 

I 2abs(^,«) J*,A dA, (Al) 

'o 

and an energy release (cooling) due to emission given by 

47ra^ (2abs(d,fl)?r6^(Pd)dd, (A2) 

where J^, J is the mean intensity of the local radiation field, 6 ^ is the 
Planck function. We could easily obtain (dEldt)^],^ = {dEjdt)^^ 
if the surrounding medium is optically thin in the wavelength 
range where the energy/radiation is absorbed/released and colli- 
sional heating is negligible. The optical depth of the medium af¬ 
fects the energy absorption but has little effect on the re-emission 
at long wavelengths, i.e., where emission can be described using 
the Rayleigh limit. Thus, in that limit and in thermal equilibrium 
we have 





0abs(d, d) JdA 



Q:bM)BA(Ti)dA, 


(A3) 


where = (2ab.s/fl- The emission at long wavelengths tends to 
follow a power-law in A. Replacing with Q'^ (/i//lo)“^ we may 
then write 


3*(a) ^ (Q,bs)* J* ^ BAT,)dA, (A4) 

where we have defined the mean value 

(Gabs)* = — r Q-Jbsi4,a) Ji,jdA, = C J^, ,ldA. (A5) 
7* Jo Jo 

7* is the wavelength-integrated mean intensity (regarded as con¬ 
stant in the present context) and f} € [ 1 , 2 ] for all common types of 
dust. The constant Go can be determined by considering Gabs^'^o^ 
for a specific dust type which conforms to a power-law approxima¬ 
tion with a constant fi. Now, the integral over wavelength on the 
right-hand side can be evaluated analytically, which yields 

^ =«4;rGo4<TsBrr^ (A 6 ) 

where cTsb is not the usual Stefan-Boltzmann’s constant, but the 
corresponding number for a modified blackbody. In case /3 = 0 
we would have (Tsb equal to the usual Stefan-Boltzmann constant. 
By Eq. 1 IA 6 I 1 it is clear that grain temperature should be weakly 
(anti-)correlated with grain size, though this correlation depends 
slightly (but not much) on the heating source and the validity of the 


appro ximations used above of course (cf. Table 3 in lPraine & Lee! 
I l984h . Note that an equilibrium model such as the one above is 
strictly valid only locally. 

If heating is due to short-wavelength radiation, we are close to 
the grey-absorption limit (particles are large compared to the wave¬ 
length) in which case (Gabs)* = 1 and thus 3* - J* = constant. 
Taking Eq. ( IA 6 b at face value, we have then a simple power-law of 
the form 

Ti{a)^TA-\ . (A7) 

\aoj 

Our aim is to arrive at a temperature distribution and for that we 
need some information abou t the GSD. A natural ansatz is th e 
canonical MRN distribution jMathis. Rumpl & Nordsiecklll977h . 
n(a) DC for which the GSD in terms of mass is if {a) oc . 
We thus have 


H'fT'd) = 


dTi 


dMi 

da 


da 






(A 8 ) 


In this special case (thermal equilibrium, MRN distribution, grey- 
absorption limit) the temperature distribution is simply also a 
power-law, which serves well as a first approximation of the func¬ 
tional form of W{T/J since 3* is typically only weakly dependent 
on the grain ra dius a. A similar derivation o f the power-law above 
can be found in lUi. Goldsmith & Xid ( Il999l) . 

In a more detailed picture the grey-absorption limit may not 
be strictly applicable and the slope of the temperature distribution 
is therefore likely steeper than in Eq. ( IA8b and may also deviate 
from the simple power-law form above, i.e., Tq would in such case 
be a function of a. However, regardless of whether 3* = 7* is a 
constant or not, Eq. l IMl l tells us that we must have ^ » 0 as 

dTi/da —> 0. That is, absence of a temperature-size relation is 
not compatible with a the existence of grains having a range of 
sizes and temperatures, unless both heating and cooling take place 
at long-wavelengths, in which case the grain radius a is cancelled 
out in Eq. IA3l l. 


APPENDIX B: COLLISIONAL HEATING 

In case heating is d ue mainly to collisions with ambient gas p arti- 
cles, we may write jPwek & WemeJl98ll : lDwekll 198611 19871) 

ldE\ I32 kV^^ , 

(■*) 'AdTA, (Bl) 

where k is Boltzmann’s constant, is the electron mass, the 
electron density, the electron temperature and h{a, Te) is a unit¬ 
less function describing the efficiency of energy d eposition. By def¬ 
initio n, h = I when the efficiency is maximal tsee lPwek & Werne3 
Il98lh . The fact that the heating rate is inversely proportional to 
the particle mass in the expression above, explains why collisions 
with electrons should be more important than collisions with any 
other gas particle, because the electron mass is very small in com¬ 
parison with, e.g., the proton mass. Since the cooling rate is the 
same as in the case of radiative heating, i.e., cooling is still due to 
long-wavelength radiation, and adopting once again a power-law 
approximation for G^bs I*’*'? wavelengths, we obtain 

1 «e h{a, Te) === a ^ (B2) 

where all quantities are as previously defined. With /i » 1 (efficient 
energy deposition), we recover a power-law of the same form as in 
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the case of radiative heating. Thus, a power-law temperature distri¬ 
bution is a reasonable first approximation not only in case the dust 
heating is due to short-wavelength radiation, but also when it is 
due to efficient collisional heating. Another important aspect of the 
above is that it is difficult to have a dust population where all grains 
have very similar temperatures unless they also have very similar 
sizes, regardless of whether the heating is radiative or collisional. 
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